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Abstract 

We consider a spectrum-sharing cognitive radio system with antenna selection applied at the 
secondary transmitter (ST). Based on the extreme value theory, we deduce a simple and accurate 
expression for the asymptotic distribution of the signal to interference plus noise ratio at the secondary 
receiver. Using this result, the asymptotic mean capacity and the outage capacity for the secondary 
user (SU) are derived. The obtained asymptotic capacities approach the exact results as the number 
of transmit antennas N increases. Results indicate that the rate of the SU scales as log(iV) when the 
transmit power of the ST is limited by the maximum allowable interference level, while rate scales as 
log(log(iV)) if ST is limited by the maximum transmit power. 

Index Terms 

Cognitive radio, extreme value theory, mean capacity, outage capacity, rate scaling law, spectrum¬ 
sharing, transmit antenna selection. 


I. Introduction 

In order to utilize the radio spectrum more effectively, the secondary users (SUs) may be 
allowed to share the spectrum with the primary users (PUs). However, the SUs need to control 
the transmit power in order not to cause harmful interference to the PUs [H|. While research 
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literature on this topic is rich, we recall in the following works that are most relevant for this 
study. Results in (21 have shown that, under a received-power constraint, the channel capacity 
of the SU can be significantly different from the one in a conventional AWGN channel. In 
|[3ll , authors studied the ergodic capacity and the outage capacity of the SU assuming fading 
channels, where the interference from the PU to the SU was ignored. In 0, authors studied a 
single-antenna spectrum-sharing system by taking into consideration the interference from the 
PU at the secondary receiver (SR). 

The multiple antenna techniques provide an efficient tool to improve the link level capacity. 
Among others, the transmit antenna selection (TAS) scheme is attractive, since it can reduce the 
hardware costs and complexity 0. In @, authors studied the ergodic capacity of the SU using 
TAS and maximum ratio combining (MRC) at the receiver. However, the PU interference was 
ignored, and no closed-form expression was given for the SU ergodic capacity in the presence 
of the peak transmit power and interference power constraints. Authors in 0 derived the exact 
expressions for the ergodic capacity of MIMO cognitive system in terms of Fox’s H-function. 
This result, although of theoretical interests, involves complicated contour integral and is thus 
difficult to use in practice. 

In what follows, we consider a spectrum-sharing system with multiple antennas at the sec¬ 
ondary transmitter. By assuming the interference from the PU to the SU and the transmit power 
constraints at the ST, we derive the asymptotic mean capacity and outage capacity for the link 
between ST and SR by applying the extreme value theory (EVT) 0. The derived results are 
simple but accurate, and they become exact when the number of antennas approaches to infinity. 
Using these results, we obtain the rate scaling laws for the spectrum-sharing systems under the 
interference power and the transmit power limitations. 

The remainder of this paper is organized as follows. In Section HU we provide the system 
model and introduce the probability density function (PDF) and cumulative distribution function 
(CDF) of the SINR at the secondary receiver. In Section HIU we briefly review the EVT, and show 
that the SU SINR of the proposed system lies on the domain of the attraction of the Gumbel 
distribution. Then we investigate the SU mean capacity, the SU outage capacity, and the scaling 
issues in Section [TV] Section [V] presents the simulation results. The last Section concludes this 


paper. 
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Fig. 1. Spectrum-sharing system with antenna selection at the secondary transmitter. Single antenna is equipped at the secondary 
receiver, the primary transmitter and receiver. 


II. System and Channel Models 

We consider a spectrum-sharing cognitive radio system depicted in Fig. [U This model is 
similar to the one in 0, where the ST was equipped with single antenna. In our model, the 
secondary receiver (SR) is equipped with one antenna, and the secondary transmitter (ST) is 
equipped with N antennas assuming TAS technique adopted at the ST. The primary transmitter 
(PT) and receiver (PR) are assumed to be equipped with one antenna. Let g { , h, (i — 1, • • • , N), 
and q denote the channel power gains between the i-th antenna of the ST and the SR, the z-th 
antenna of the ST and the PR, and the PT and the SR, respectively. 

When the i-th antenna is used, the signal received at the SR, can be written as 

Vi = V p s9i s + a /P^qz + v, (1) 

where P s and P p are the transmit power of the SU and PU, respectively, s and z denotes 
the transmit signals of the SU and PU with unit average power, and v represents the additive 
white Gaussian noise with power a 2 . We assume that the channels experience independent block 
Rayleigh fading. Namely, the channels keep constant during each block, and are uncorrelated 
between blocks ®. Furthermore, the channels between the ST and the SR are identical and 
independent distributed (i.i.d.), and the channel between the ST and the PR is independent from 
the secondary channels. For a Rayleigh fading channel, the channel power gain x, x e { g t , h,, q }, 
follows exponential distribution with mean value x. The probability density function (PDF) of 
x is given by 

fix) = =e~ x/x , x > 0 
x 


(2) 
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and the cumulative distribution function (CDF) is 

F(x) = 1 - e~ x/m , x>0. 


( 3 ) 


The selected antenna element of the ST provides the maximum output SINR at the SR. The 
SR then sends back the index to the ST through an error-free feedback channel. In order to limit 
the interference caused by the SU, an interference power constraint, denoted by Q, is adopted, 
which limits the instantaneous interference generated by the SU to the PU. In addition, P max 
denotes the maximum transmit power of the ST due to, for instance, hardware limitations or 
regulations. On limited power, now, the transmit power of the ST using the Pth antenna can be 
represented as 0], lilOl 

P® =min j^,P max }. (4) 

Accordingly, when the Pth antenna of the ST is used, the received SINR at the SR is represented 
as 


7 i = mm 


Qdi 


P„ 


hi ( P p q + a 2 )' P p q + a 2 

The CDF of 7 * has been deduced in 0 eq. (15)] as 


P 7i (P) = 1 


Qhi 


Qai 


-=e 


PpQ xPpqh^ 


xP p qhi 


( 5 ) 


x r 


0 , 


X 


1 — e 


i/jP]n;lx 

Q 


+ 


P P Q 

9iP 1 


xh, 


( 6 ) 


xP p q + g t P xr 


- g 9 i Pmax 


where T(a,6) = J,°° e _t f a_1 df denotes the incomplete Gamma function ifTTl eq. 8.350-2]. If the 
TAS technique is used, then the antenna of the ST providing the maximum SINR at the SR is 
selected such that 


7max = max 7;, i = 1, • ■ ■ , N, (7) 

l 

where 7 * is given by ©. The CDF of the received SINR at the SR can be written, using order 
theory [|T 2 ll . as 

P^(x) = [F^(x)} N . ( 8 ) 


To the best of our knowledge, there is no explicit expressions for mean and outage capacities 
when using the CDF ([ 8 ]) due to its algebraic complexity. In order to obtain tractable analysis 
and gain insight for the considered cognitive radio system, we deduce asymptotic expressions 
for the mean capacity and outage capacity using extreme value theory. 
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III. Implications of the Extreme Value Theory 


Let X 1 , ■ ■ ■ ,X N denote the i.i.d. random variables, and let F x {x) be the underlying cumulative 
distribution function. We denote by M x the maximum of {2Q} 1<i<iV . If there exist real constants 
a N and b N > 0, and a non-degenerate limit distribution G(x) such that 


M x — a x d , 
hm --->■ G(x), 

N^oo b N 


(9) 


then G(x) is one of the extreme value distributions: Gumbel, Frechet, or reverse-Weibull |[T3l . 
We say that F x (x) belongs to the maximum domain of attraction (MDA) of G{x). 

Before using the extreme value distributions to approximate F 7max ( x) in ([8]), we recall a lemma 
proved in [HI, [fl4l which gives a sufficient condition for a distribution belonging to the MDA 
of the Gumbel distribution. 

Lemma 1: [[8] Theorem 2.7.2] and [14, Lemma 1]. Let u{F x ) = sup {x : F x < 1}. Assume 
that there is a real number xi, such that for xi < x < co(F x ), f x (x) = F' x {x) ^ 0 and F'i(x) 
exists. If 

d n - Fv(x\ 

= 0 , 


lim — 

x- h*i(Fx) ax 


fx(x) 


then F x {x) lies on the MDA of the Gumbel distribution ([9]) with the CDF expressed as 

G(x) = e" 


x — apj 
bpj 


( 10 ) 


( 11 ) 


where the normalizing constants a x and b x are determined by 


a x = F 


-l 

x 


bn = F 


-l 

X 


1 

1 — 

N 

1 

Ne 
x X ( 


( 12 ) 


— a N . 


Here e is the base of the natural logarithm, and F x (•) denotes the inverse of the CDF. 

Using Lemma |T] we prove that the CDF of SINR TA (•) lies on the MDA of the Gumbel 
distribution. 

Proposition 1: Assume that the transmit power of the secondary transmitter is given by ©. 
Then the SINR distribution F lt (■) of the SR with transmit antenna selection belongs to the MDA 
of the Gumbel distribution. 

Proof: see Appendix [A[ ■ 

Using Proposition [0 we obtain approximations for the mean capacity and outage capacity of 


the SU. 
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IV. Asymptotic Performance Analysis 


A. SU Mean Capacity 


The instantaneous rate for the SU can be represented as 


R, 


-max 


max R ,, 

l<i<N 


(13) 


where F, = log (1 + 7 ^), \/i = 1, • • • , N is the instantaneous rate using the i- th antenna. Accord¬ 
ing to the limiting throughput distribution (LTD) theorem proved in Ifl4l . the distribution of R, be¬ 



longs to the MDA of the Gumbel distribution. Consequently, we have lim/v-xx> 
and G{x) is expressed in (ITTI) with the following normalizing constants 



(14) 



where F~. 1 (•) denotes the inverse function of P 7 . (■) in ©. Based on the Lemma 2 in lfl4ll . we 
find that the expectation of converges, and the mean capacity of the SU C s = E [f?, max ] 

can be thus approximated by integrating F max over (HU) as 


C s ~ + b N E 0 , 


(15) 


where E 0 = 0.5772 ... denotes the Euler constant. According to Proposition |T| the approximation 
is tight in the asymptotic regime N — > oo. 

Based on (fl5l) . we study the behavior of the SU mean capacity in two limiting regimes. First, 
we consider that the maximum transmit power constraint approaches infinity (P max — > oo). And 
the second case is that the interference power constraint is much larger than P max , i.e., Q F> P max - 
We name these two scenarios as interference power limited regime (IPLR) and transmit power 
limited regime (TPLR), respectively. 

1) Interference power limited regime: When P max — > oo and the z-th antenna is selected, we 
rewrite the asymptotic CDF © of the SINR at the SR as 



where c p = e! 7 ' 1 ^ Pvq )\ ^0, -jf=J /P p q represents the effects of the additive noise and the PU 
interference. The inverse CDF of the SINR is given by 



( 17 ) 
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Using (fT4l) and (IT7l) . the normalizing constants of the Gumbel distribution G(x) are obtained as 

Qg 


IPLR 

l N 


and 


= log ( 1 + c p ^-N ) , 


ip, T, i 1 + Cp^-Ne 

b 1RLR = log -- 

N 81 1 + c pfN 


(18) 


(19) 


where b^ LR —* 1 as N —* oo. Thus, the asymptotic mean capacity (IT5l) becomes 

Qg 


C = log ( 1 + c p -=-N ) + E 0 , 


( 20 ) 


which indicates that the SU mean capacity scales as log(iV) in the interference power limited 
regime. 

2) Transmit power limited regime: When 0 P max and the z-th antenna is selected, we 
rewrite the asymptotic CDF © of the SINR at the SR as 


R 


7i,TPLR 


x) = 1 - 


Cq Pmax 9 

xP p q + P max g 


g Pmax 9 ^ 


( 21 ) 


where c q = 1 — e p// ' and p = Q/P max . The inverse CDF of the SINR can be obtained as 


p- 1 

- r 7i ,TPLR 


r 1 Pmax9 . , 

[: v) ~ ——W 


c q e p p q 


<y 


Pm ax.9 

(i - y)P P q/v 2 1 P P q 


( 22 ) 


where W(x) denotes the Lambert W-Function satisfying W(x)e yv ^ = x. Thus, using (fl4l) and 
|), the normalizing constants of the Gumbel distribution G(x) are obtained as 


< = log 1 + 


( P Cq C Hp q 


-w 


P p q/c 


b'^ K = log I 1+ I eNCq ^ 


P p q/<J 2 


PuiAyfJ 

P p q 

Pmax.9 

PpQ 


(23) 


TPLR 
2 N 1 


(24) 


where //y LR = O (1 /N) as N —* oo, and 0 represents a term of order ( z). Therefore, the 
SU asymptotic mean capacity ([131) becomes 


^yTPLR 


= a 


N 


+ b 


N 


(25) 


which admits an explicit expression and avoids the operation of inverting the SINR CDF function 

©. 














When the interference from the PU can be ignored, i.e. P p q —* 0, (I2TI) becomes 

xcr 2 

P^TPLRW ~ 1 - C q e 
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(26) 


which is the approximation of the CDF of the SINR at low interference regime. Then we obtain 
the inverse CDF as 

(27) 

Substituting (l27l) into (fT4l) . we have the normalizing constants 


a 


TPLR 

iV,low 


= log 1 + 


PnuxfJ 


<J 


lOg(CgA^) , 


1 

77 / ^2 


f TPLR 

-'A 1 ', low 


= log 1 + 


PmaxS / ® 


1 + log (c q N) 


Substituting (l28l) and (|29|) into (fl5l) . the asymptotic mean capacity is obtained as 


(28) 

(29) 


^TPLR 

°low 


TPLR , /TPLR p 
u iV,low "r "TV.low^O) 


(30) 


where = O (1 /N) as N —> oo. This shows that the SU mean capacity scales as loglog(iV) 
in the transmit power limited regime with low PU interference. Note that in the TPLR with low 
PU interference, the considered model £0) reduces to the conventional point-to-point communi¬ 
cation. In this case, the derived result (l30l) is in line with [15] eq. (38)]. 

Remarks: The results suggest that the TAS provides different power gains in these two 
scenarios. In the TPLR, (l30l) shows that increasing the number of transmit antennas may not be 
an efficient way to achieve higher mean capacity, while it is effective in the IPLR as shown in 
£ 20 ]). 


B. SU Outage Capacity 

Given the rate r (nats/s/Hz), the SU outage probability after transmit antenna selection can 
be defined as 

^out(r) = Pr {-R max < r} , (31) 


where R m;iX is given by (fl3l) . The outage capacity can be interpreted as the maximum rate such 
that F out < e, where e is the defined outage threshold. For monotonically increasing function 
F out (r), we have 


Gout — 


^(e). 


(32) 








According to the LTD theorem [fl4l1 . the asymptotic distribution of R max can be approximated 
as CD- Thereafter, we have the following proposition for the SU outage capacity. 

Proposition 2: With TAS, the asymptotic outage capacity of the SU is given by 


Cout « a N -b N log log 



(33) 


\e j 

where the normalizing constants a ; y and b N are the normalizing coefficients given in (fT4l) . 

Proof: The proof of this proposition is straightforward by substituting (fTTl) into (l32l) . ■ 

Based on Proposition [2] we investigate the SU outage capacity scaling properties in case of 
IPLR and TPLR, respectively. 

1) Interference power limited regime: When P max —> oo, the normalizing constants a IRLR and 
b% LR are computed in (fl8l) and (fl9l) . respectively. Following (l33l) . the asymptotic SU outage 
capacity in IPLR yields 


^out = log ( 1 + c p~=~N ) ~ log 


Os ; 

h 


'l +Cp QjlNe\ 


1 C' 


Ul__ 

*¥ N ) 


log log 


(34) 


Recall that 6T LR —> 1 when A r ^ oc, we have 


log ^1 + c P^ N ^j ~ lo £ lo g Q (35) 

= C lPLR - E 0 - log log , (36) 

where 6' IPI R is given in (l20l) . We can see that the SU outage capacity scales in the same 
manner as the SU mean capacity. When given the outage threshold e, there exists a constant gap 
E 0 + log log (i) between C' ( IRRR and C IPLR . 

2) Transmit power limited regime: When 0 P Iriax , following (l33l) the asymptotic SU outage 
capacity in TPLR yields 

C™ = «™-f-» LE loglog(i), (37) 

where the normalizing constants a^ LR and /;(f LR are computed in (1231) and (1241) . respectively. 
We can see that the SU outage capacity scales the same way as the SU mean capacity, i.e. 
C™ ~ log log (. N) as —> cxo. It is also shown that the asymptotic SU outage capacity does 
not depend on the outage threshold e. However, the convergence speed is pretty slow due to 

C*™ ~ l°g l°g (N). 
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When the interference from the PU can be ignored in TPLR, i.e. P p q —> 0, the asymptotic 
SU outage capacity ([371) yields 



(38) 


where a™ R and are given by (1281) and (|29|) . respectively. 

Remarks: The number of antennas has different influence on the asymptotic outage capacity 
of the SU. In the TPLR case, i.e. P max -C Q, the SU outage capacity increases logarithmically 
with N and the gap between the SU mean capacity and outage capacity vanishes. However, this 
is not for the IPLR scenario, i.e. P max oo, that the outage capacity has a much slower scaling 
rate, log log(TV). In addition, there is always a gap between the outage capacity C^ R and mean 
capacity U IPLR . 


V. Simulation Results 


We present the asymptotic results with respect to N for the SU mean capacity and the outage 
capacity, assuming that the ST employs transmit antenna selection and the transmit power is 
restricted by the peak interference power constraint and the maximum transmit power constraint. 
The rate scaling behaviors of the mean capacity and outage capacity are also discussed. The 
mean channel power gains are set to 1. The noise power per Hz at the SR is set to —10 dB, 
i.e., a 2 = 0.1, where we assume that the maximum average received signal to noise ratio at the 
SR is greater or equal to 1, i.e., E [P max gi\ /a 2 > 1. Let INR = P p /a 2 denote the interference 
to the noise ratio, SIRq = Q/P p and SIRp = P p be the transmit signal to the interference 
power ratio, and SNRq = Q/cr 2 and SNR P = P max /cx 2 be the transmit signal to the noise power 
ratio. In addition, PQR = P max /Q is the maximum power to the interference constraint ratio. 

Fig. |2 shows the SU mean capacity versus the interference power constraint 0 when the 
number of the ST antennas is N = 4, 10, and 20. The PU transmit power is set to P p = 0 dB, 
and the maximum transmit power P max of the SU is 0 dB. We also plotted the asymptotic mean 
capacity of the SU according to (fL5l) . and compared it to the simulation results generated by 10 6 
channel realizations. As the interference power constraint Q increases, the SU mean capacity is 
improved until the maximum transmit power constraint P max becomes dominating. We can see 
that the asymptotic results using EVT are reasonably accurate, which validates our approach to 
characterize the SU mean capacity using simple and explicit expressions. 
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Fig. 2. SU mean capacity versus the peak interference power constraint Q with P max = 0 dB and P p = 0 dB. 

Fig. [3] shows the the SU mean capacity as a function of the number of ST antennas in 
the interference power and the transmit power limited regimes, respectively. In Fig. [3a] the 
maximum transmit power constraint is set to P max = 30 dB3> Q. In this case, we calculate 
the SU mean capacity using IPLR approximation, namely (fl5l) with the normalizing constant 
a IRLR and 6 IRLR given by (fl8l) and (fT9l) . respectively. In addition, we show the SU rate scaling 
(l20l) in IPLR. Numerical results show that the approximation and scaling yield good agreement 
with simulations, especially when PQR = 35 dB and SNRq = 5 dB. Furthermore, the scaling 
behavior of the SU capacity is shown to be log(iV) as predicted by (l20l) . In Fig. [3b] we assume 
that PQR = {—10, —20} dB. The SU mean capacity is computed using TPLR approximation 
(1251) . and the SU rate scaling is obtained via (l30l) . In TPLR with low PU interference, e.g. 
INR = —10 dB, the SU capacity scales as log (log (IV)). Results show that in TPLR, it is more 
effective to improve the SU capacity by increasing the transmit power than by using more 
transmit antennas. However, in IPLR if the SU is subject to a lower interference limit, the mean 
capacity can be maintained by using more transmit antennas. 

The SU outage capacity is plotted as a function of the interference constraint Q in Fig. [4] 









12 



(a) IPLR 



(b) TPLR 


Fig. 3. SU mean capacity versus the number of ST antennas N. Lines with circles denote INR = —10 dB; lines with diamonds 
denote INR = 10 dB. (a) IPLR: PQR = {35,40} dB. (b) TPLR: PQR = {-10, -20} dB. 























13 



Fig. 4. SU outage capacity. P max = 0 dB, P p = OdB, and e = 10%. 

where the outage probability is set to 10%. Although the EVT approximation is not very accurate 
for N = 4 ST antennas, it gives an estimate trend of the SU outage capacity. As the number of 
antennas increases, the accuracy of the EVT approximation improves. In addition, for N = 20 
the approximation for a large range interference constraint, Q > — lOdB, is acceptably accurate 
for understanding such systems. 

In order to gain some insight into the asymptotic behavior of the SU outage capacity, we depict 
in Fig. |5| the SU outage capacity in the transmit power limited and interference power limited 
regimes, respectively. The asymptotic SU mean capacity is also plotted for a comparison. As the 
number of ST antennas increases, the asymptotic outage capacities for various cases are accurate 
compared to the numerical simulations. The scaling properties of the SU outage capacity are 
validated. In Fig. |5al the SU outage capacity in IPLR behaves as log(iV) when the number of 
ST antennas N approaches to infinity. In this case, there is a constant gap E 0 + log log (1/e) 
between the outage and the mean capacities as predicted by (l36l) . Yet, the SU outage capacity 
in TPLR behaves as loglog(lV) as shown in Fig. |5b] which can be proved using (l38l) . 
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Fig. 5. SU outage capacity versus the number of ST antennas N. Lines with circles denote INR = —10 dB; lines with diamonds 
denote INR = 10 dB. The peak interference power constraint Q = 5 dB. Outage threshold e = 10%. (a) IPLR: PQR = 25 dB. 
(b) TPLR: PQR = {-5, -15} dB. 
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VI. Conclusion 

In this paper, we applied the extreme value theory while studying the performance of a 
spectrum sharing system in terms of the SU mean capacity, the SU outage capacity, and the rate 
scaling characteristics. The SU mean capacity and outage capacity have two scaling behaviors 
with the number of available transmit antennas: in the interference power limited regime (IPLR), 
capacities scale as log(iV), while they scale as log(log(iV)) in transmit power limited regime 
(TPLR). These results indicate that the transmit antenna selection technique provides different 
power gain in considered two scenarios. The obtained accurate approximations provide us a 
better understanding on the performance of the spectrum-sharing systems with multiple transmit 
antennas. 


Appendix 

Derivation of Proposition HI 

The associated PDF becomes, through taking the derivation to (x) with respect to x and 


expanding at x —> oo, 


fji( x ) — e 


o 


Q 


\ — g -Pmaxfr' 


P p qx 


+ 0 


x 


(39) 


where O (z) m represents a term of order (z)’". The expansion at x -x oo of the secondary 
derivative, f (x), can be expressed as 


<7 


f'(x) = e ^ 


e r "»x“ i / -| 

+ 0 - 

U| i;i \/V / (//•/' \ ■!' 


It it obvious, from (l39l) and (l4fl . that lim^oo / 7 . (x) = 0 and Hindoo [—f { 
respectively. We need to show, after taking the derivative in (fTOl) . that 

■(1 -F 7i (t))/4(x 


x 


lim 

x—>oo 




= -1 


(40) 
= 0 , 

(41) 


( x ) is the PDF of 7 j, thus we have 


lim [F 7 . ( x ) — 1] = 0 
£—>•00 1 


( 42 ) 
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Also, we have lim^oo [/ 7i (x)] = 0. Therefore we obtain, according to L’Hospital’s rule, 

(1 -A.M)' 


lim 

#—>•00 




Thus, link, 


fn( x ) 

' U(x) 

x ^°° L -f^( x ) 

l-F^Ax) 


= lim 

x —>-oo 


/ 4 ( x ) 


= lim 


fa O) 


7* 

or lirrij^oo 


= lim ( + o(- 

x—roo \ (J 2 \ X 


PmaxQ 


(43) 




f~n (*) 


lim 

x—>oo 


(l- F^x^f'Ax) 


both have a finite limit. Therefore, we have 


= lim 

a?->-oo 


w 

1 ~ F li^) 

fyX x ) 


r i 

lim TtT 
x ^°° v fa \ x ) 


= -1 


(44) 


This completes the proof. 
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